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' ABSTRACT 

Formulas relating to the elastic behavior of axisym- 
metrically loaded cylindrical and spherical pressure vessels 
have been collected from various sources. Their presenta- 
tion is unified by the employment of a uniform system of 
notation. Loadings include interior and exterior pressuriza- 
tion, radial and axial temperature gradients, axial loads, 
and centrifugal force fields. Some cases of buckling are 
treated, and shrink assembled multilayer vessels are 
included. A large bibliography is listed. 
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NOTATIONS 



(Dimensions are indicated in brackets) 

In the formulas which appear in this thesis, great care 

is exercised to assure that superscripts and exponents can 

2 

not be confused. Thus, we write n x n = (n) rather than 

p 

n , since the latter could be interpreted as the value of 
n which pertains to section two. However, the expression 
1 - (v) occurs so frequently that we write 1 - v , 
omitting the parentheses for simplicity. 

ENGLISH LETTERS 



b [L] 

B [FL -4 ] 
C [-] 

c [-] 

D [FL] 

D c [FL] 

D g [FL] 

D p [FL] 

d [L] 

d c CL] 



Distance from the end of cylindrical pressure 
vessels to any selected point. 

A parameter used only in Section XLIII. 

2 

Parameter = cosh(B x). C = cosh(8 5-). 

c c 

2 

Parameter = cos(8 s). c = cos(B A). 

Shell flexural rigidity in general. 

Flexural rigidity of cylindrical pressure 
vessels _ E(t c )^ 

' 12(l-v 0 2 ) 

Flexural rigidity of spherical pressure 
vesseis ^>3 

12(l-v 2 ) 

E(t)3 

Flexural rigidity of plates = E__ — . 

12(l-v p 2 ) 

Diameter of pressure vessels in general. 

Mean diameter of cylindrical pressure vessels. 
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ci 



CO 



si 



so 



E 

E, 



cr 






i 

M 



[L] 

CL] 

CL] 

[L] 

CL] 

[FL -2 ] 

[FL -2 ] 



[fl -2 ] 

[F] 

[F] 

[LT -2 ] 

[FL -1 ] 

[FL -1 ] 



C-] 



[-] 



[FL -3 ] 



[FL -3 ] 



[L] 

[F] 



Inner diameter of cylindrical pressure vessels 
Outer diameter of cylindrical pressure vessels 
Mean diameter of spherical pressure vessels. 
Inner diameter of spherical pressure vessels. 
Outer diameter of spherical pressure vessels. 
Modulus of elasticity. 

Modulus of elasticity of i^* 1 cylinder of 
multilayer pressure vessels. 

Tangent modulus of elasticity. 

Axial tensile load on pressure vessels. 

Critical value for axial buckling force. 

Gravitational acceleration. 

Edge Shear force (per unit length) in general. 

Edge shear force on pressure vessels in which 
the subscript j denotes the junction number 
where the edge shear force exists. 



Stress index in general = 



(El) * 

^2t ; 



Stress index in which subscript x denotes the 
type of pressure vessel; superscript Y denotes 
the type of loading, j denotes the junction 
number where the stress is evaluated. 

Parameter for cylindrical pressure vessels 
_ Et 



<V 



Parameter for spherical pressure vessels 
Et s 

CR S ) 2 

Length of cylindrical pressure vessels. 

Edge bending moment (per unit length) in 
general. 
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M. 



M 



cr 




R 

R 

c 

R 

R 

c 

R 



ci 



co 



si 



co 

4 



r 



[F] 

[FL] 



C-] 

[FIT 2 ] 

[FIT 2 ] 

[FL -2 ] 

[FL" 2 ] 

[FL -2 ] 



[FL" 1 ] 

[L] 

[L] 

[L] 

[L] 

[L] 

[L] 

[L] 

[L] 



CL] 



Edge bending moment on pressure vessels in 
which subscript j denotes the junction number 
where the edge bending moment exists. 

Critical value for torsional buckling moment 
on cylinder. (Note dimensions are [FL] 
rather than [F].) 




Pressure in general. 

Critical value for external pressure. 

Internal pressure. 

External pressure. 

Successive interface pressures of pressurized 
multilayer vessels with shrink-fit, in which 
subscript x denotes the type of pressure 
vessel, k denotes inner or outer surface, and 
superscript i denotes the layer number. 

Axial tensile force per unit circumference 
F 

2irR 

c 

Radius of pressure vessels in general. 

Mean radius of cylindrical pressure vessels. 

Inner radius of cylindrical pressure vessels. 

Outer radius of cylindrical pressure vessels. 

Mean radius of spherical pressure vessels. 

Inner radius of spherical pressure vessels. 

Outer radius of spherical pressure vessels. 

Radius of each pressure vessel member of 
multilayer pressure vessel in which subscript 
x denotes the type of pressure vessel, k 
denotes inner or outer radius, and i denotes 
the layer number. 

Distance in radial direction (subscript r 
indicates radial direction) 
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s 


[-] 


Parameter 


= sinh(B x). : 

0 


= sinh(3 c &) 


s t 


1 — 1 
C\J 

1 

pH 
1 1 


Equipment 


allowable tensile 


stress . 


s 


C-] 


Parameter 


2 

= sin(B x). s 


sin(£ c Jl) . 



In the immediately following definitions "temperature" 
means increase in temperature over a uniform initial 
temperature. 




T 

T i 

T o 

t 



t 



c 



t 



e 



t 



cn 



t 



P 



t 

s 



t 



t 



u 



u 



i 



[°F] Temperature, as a function of radius. 

• [°F] Temperature difference in radial direction, 
positive if the outside is hotter than the 
inside . 

[°F] Average temperature through wall of vessel 
as a function of axial position. 

[°F] Temperature of inner surface. 

[°F] Temperature of outer surface, also axial 
difference in temperature in Section XLV. 

[L] Thickness of pressure vessels in general. 

(Subscript indicates tangential or 
circumferential direction, except see t. 
below. ) 

[L] Thickness of cylindrical pressure vessels. 

[L] Thickness of ellipsoidal pressure vessels. 

[L] Thickness of conical pressure vessels. 

[L] Thickness of plate. 

[L] Thickness of spherical pressure vessels. 

[L] Thickness of toroidal pressure vessels. 

[L] Radially outward deformation of pressure 

vessels. 

[L] Deformation of pressure vessels in which 

subscript i denotes pressure vessel number 
to be taken. 
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X 


[L] 


Distance in axial direction. (Subscript x 
indicates axial direction.) 


y 


[L] 


Distance in y direction. 

(Subscript y indicates y direction.) 



GREEK LETTERS 



a [(temperature)”' 1 '] 

Coefficient of thermal expansion. 



*0 


1 

1 1 


6 s 


CL" 


Y 


1 — I 

t -1 


6 


[L] 


6 * 


1 — I 
1 1 


6 o 


1 — 1 
1 1 


{ ci 


[L] 


6 co 


1 1 

1 1 




[L] 


4 si 


[L] 




[L] 


S 1 

X 


1 1 

1 1 



Parameter for cy lindrical pressure vessels 
Et_ ,'|3(l-v/) 



‘•W 



<V 2 (t e )2 



Parameter for sp herical p ressure vessels 



Et 



"W 



1 3(3 -V> 

(R s ) 2 (t s ) 2 



Specific weight, (i.e., weight per unit volume). 
Deflection or displacement in general. 



Radial interference (as manufactured) prior 
to assembly of layers of multilayer pressure 
vessels . 



Radial deflection or displacement of cylindrical 
pressure vessels. 

Radial deflection or displacement of inner 
surface of cylindrical pressure vessels. 

Radial deflection or displacement of outer 
surface of cylindrical pressure vessels. 

Radial deflection or displacement of spherical 
pressure vessels. 

Radial deflection or displacement of inner 
surface of spherical pressure vessels. 

Radial deflection or displacement of outer 
surface of spherical pressure vessels. 

Radial deflection at the edge (end) of pressure 
vessels in which subscript x denotes the type 
of pressure vessel, superscript i denotes the 
junction number where edge deflection occurs. 
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6 x Yj C * ] 



[-] 



The edge deflection influence coefficient of 
pressure vessels in which the subscript x 
denotes the type of pressure vessel, super- 
script i denotes the junction number where the 
edge deflection coefficient is to be taken, 

Y denotes the type of loading, j denotes the 
junction number where the loading exists. 

[*] indicates that dimension depends on 
the dimension of the corresponding (Y) 
loading. 

Unit strain. 



6 [— ] 



C— ] 



efj [•] 



\ 



e 



C-] 



[-] 



[-] 



v [-] 

v ± C-] 



Rotation of edge of pressure vessels in 
general. (All angles are measured in radians 
except when the symbol ° is added; e.g., 0°.) 

Edge rotation of pressure vessels in which 
subscript x denotes the type of pressure 
vessel, superscript i denotes the junction 
number where the edge rotation is to be taken. 



Edge rotation influence coefficient of pressure 
vessels in which the subscript x denotes the 
type of pressure vessel, superscript i denotes 
the junction number where the edge rotation 
coefficient is to be taken, Y denotes the 
type of loading, j denotes the junction 
number where the loading exists. [*] indicates 
that the dimension depends on (Y) loading. 



Parameter for cylindrical pressure vessels 

■ iW 2 ) <^> 2 

Parameter for spherical pressure vessels 

■ ^ U - v s 2 ) (^) 2 

s 



Parameter = 



m 



<3 t 
c c 



H, 

H, 



Poisson’s ratio. 



■f* V* 

Poisson’s ratio of 1 cylinder of multilayer 
pressure vessels. 



5 [-] Parameter in Section XIX. 
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p 


[FT 2 L" 4 ] 


v 

Mass per unit volume = — . 

g 


a 


1 1 

OJ 

1 

1 1 


Stress In general. 


a cv 


1 — 1 
OJ 

1 

1 1 


Critical value for compressive axial 
buckling stress of pressure vessels. 


a l 


1 — 1 
C\J 

1 

i-q 

t i 


Meridional stress. 


a 

r 


1 — l 

C\J 

1 

I I 


Radial stress. 


o 1 

r 


I — i 

OJ 

1 

1 1 


4~ V'i 

Radial stress in i^ layer of multilayer 
vessel . 




1 — 1 
OJ 

1 

1 1 


Tangential or circumferential stress. 




[fL' 2 ] 


Tangential or circumferential stress of 
ith layer of multilayer vessel. 




1 — 1 
00 

1 

& 
1 1 


Larger hoop or circumferential tension at 
radius R^ o in multilayer pressure vessels. 


<Vs 


[PL' 2 ] 


Smaller hoop or circumferential tension at 
radius Ri Q in multilayer pressure vessels. 


a 

x 


[PL' 2 ] 


Axial stress. 


a 1 

X 


1 1 

00 

1 

1 1 


t h 

Axial stress in l layer of multilayer 

vessel . 


a 

yp 


[FL -2 ] 


Yield stress. 


T 


[FL -2 ] 


Shear stress in general. 


T 

cr 


1 1 

00 

1 

1 1 


Critical value for buckling shear stress due 
to torsional moment. 


$ 


C-] 


Angle in circumferential direction measured 
from any selected point. 


* 


C-] 


Angle in longitudinal direction measured 
from a plane perpendicular to the axis of 



pressure .vessel. 
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I . INTRODUCTION 



This thesis presents a collection of formulas applicable 
to mechanical (structural) analysis and design of axi- 
symmetric pressure vessels under axisymmetric loadings. 
Emphasis is on cylindrical vessels, but some results are 
given for spherical vessels. Furthermore, a variety of 
end closures are considered in connection with long thin- 
walled cylindrical pressure vessels. 

No claim is made for any originality of content. For 
the most part the results presented herein may be considered 
to be a part of the general or common knowledge of the 
mechanics of pressure vessels and they are given in many 
references. Generally, one or more literature citations 
are indicated herein pertaining to each set of results. 

The literature contains several instances of what can be 
considered as errors by authors and compilers of good repu- 
tations, but which we choose to consider as misunderstandings 
resulting from imperfect specification of auxiliary condi- 
tions. For this reason we have refrained from pointing out 
such cases in the literature but have attempted to assure 
that our own reporting of results is not subject to such 
misunderstandings . 

No attempt is made herein to give derivations of results 
which are dealt with extensively elsewhere except that 
Appendices 1 and 2 hereof do provide, in what we believe to 



15 



be a very compact and useful form, basic analyses for semi- 
infinite and finite (respectively) thin-walled cylinders, 
results from which are employed extensively in the body 
thereof. Also, Appendix 3, contributed by Professor Brock, 
provides, in compact form, an analysis which, though 
elementary, we have not found elsewhere. 

Also, no claim is made for completeness. The literature 
is so very extensive and is growing at such a rate that it 
would be quite impossible to report all significant and 
useful results. Time and space limitations have thus 
limited the scope of what is here presented. What does 
appear, then, reflects the interest of the author who has no 
wish to defend his choices or to rebut any charges of 
inconsistency in selection. 

However, it is indeed hoped that this compendium will 
be found useful by- those into whose hands it falls. Most, 
or at least a great many, of the commonly encountered load- 
ing situations are treated. An attempt has been made to be 
quite detailed and complete in delineating all pertinent 
constraints and conditions. (For example, the influence of 
axial strain or axial loading is explicitly considered, thus 
tending to avoid the errors or misconceptions which frequent- 
ly arise when this matter is evaded.) In some -cases ade- 
quately treated in the literature it is nevertheless diffi- 
cult to obtain a clear picture of alternate modes of behav- 
ior and the circumstance under which each may occur. In 
these cases we hope that the broader picture may be usefully 
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evident herein. (One such case is that of torsional 
buckling of thin-walled cylinders.) Finally, the literature 
citations which appear in the body hereof and the somewhat 
more extensive listings in the bibliography should prove 
useful in guiding the reader to sources of information he 
may need but which does not appear herein. 

One general remark should perhaps be appended. Through- 
out it is assumed that the material behaves in a linear 
elastic manner and that the resulting strains are sufficiently 
small that second-order effects can be neglected. 

It is also assumed, except in cases specifically dealing 
with buckling, that deformations do not significantly alter 
the loading configuration. Briefly, except for cases of 
buckling, strictly linear behavior is postulated. 
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II. MEMBRANE THEORY 



This analysis is thoroughly treated in the pressure 
vessel literature [4, 10, 17, 19, 25, 26, 27, 32,33, 3*»]. It is 
based upon considering a thin plate or sheet of revolution 
symmetrical about axis x-x (Pig. 1). 







Pig. 1 Element' of Axisymmetric membrane 

If the wall thickness t of such a shell is small compared 
with its principal radii of curvature, the wall will have 
negligible bending resistance and acts as a membrane in which 
the stress is tangential to the middle surface of the wall 
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and uniformly distributed across its thickness. These 
stresses are called membrane stresses and can be found from 
equations of statics alone. 



A basic equation describing the equilibrium of an element 
(see Fig. 1) in a direction normal to its surface is: 



Results for thin-walled cylindrical shells are reported 
by many authors [4,6,10,17,19,25,26,27,32,33,37]. It is 
assumed that the shell is long enough that the edge shear 
forces and bending moments have negligible influence on the 
section considered. 

The well-known stress and deflection equations for thin- 
walled cylindrical shell can be written as follows: 



+ s; fr-fe 



r t r t ± 



(l) 



III. LONG THIN-WALLED CYLINDRICAL SHELLS UNDER 
INTERNAL AND EXTERNAL PRESSURE 








(3a,b) 
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Formulas (a) are for the case of a general axial tensile force 
F in the material of the cylinder. If this arises solely from 
pressure, then F = (p^-p^TTF^ and formulas (b) result. 

IV. LONG THIN-WALLED CYLINDRICAL SHELLS 
UNDER INTERNAL PRESSURE 

The stress and deflection equations can be obtained 
from Equations 1,2, and 3 in Section III by putting p Q = 0. 

The* simplified formulas are 





% k 


(la,b) 




‘•fV' 


(2a, b) 


*6 2'T ft. 




(3a, b) 



V. LONG THIN-WALLED CYLINDRICAL SHELLS 
UNDER EXTERNAL PRESSURE 



Equations of stress and deflection can be obtained, as 
was done in Section IV, from equations 1, 2, and 3 in Section 
III by putting p^ = 0. The results are 






(la,b) 
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(2a, b) 




JfrK c t c 









(3a, b) 



cf. , also. Section XXXV 



VI. THIN-WALLED SPHERICAL SHELLS UNDER 
INTERNAL AND EXTERNAL PRESSURE 



Stress and deflection equations for spherical shells can 
be found In many books and articles [4,6,10,17,19,25,27,33, 
34,37]* Membrane equation 1 in Section II can be used to 
get the stress equation for thin-walled spherical shells by 
taking longitudinal and hoop radii to be equal and also 
taking longitudinal and hoop stress to be equal. Hence the 
equation can be written as follows; 




( 1 ) 



The deflection equation may be written as follows 




( 2 ) 
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VII. THIN-WALLED SPHERICAL SHELLS UNDER 



INTERNAL PRESSURE 



Stress and deflection equations can be obtained from 

equation 1 and 2 in Section VI by simply putting p Q = 0 as 

\ 

follows 




( 1 ) 




!z± < 4 , 

2 e ± i r * 



( 2 ) 



VIII. THIN-WALLED SPHERICAL SHELLS UNDER 
EXTERNAL PRESSURE 



Substituting p^ = 0 in equations (1) and (2) of Section 
VI gives 



<S'- 






>-y 

2E 



Es b 
t s r ° 


(1) 


*. ' 


(2) 



(Cf., also Section XLI.) 



IX. LONG THICK-WALLED CYLINDRICAL PRESSURE VESSELS 
UNDER INTERNAL AND EXTERNAL PRESSURE 

Thick-walled cylindrical vessels have been studied by many 
authors [4,6, 19,33 s 34,36] . When the thickness of the 
cylindrical shell is relatively large, the stress variation 
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from inner to the outer surface becomes important. (The 
analysis in this section is known by the name of G. Lame 
who first developed it.) 

The membrane theory can not be applied. If a thick 
walled cylindrical vessel with constant thickness is sub- 
jected to uniform internal and external pressure, 
deformation will be symmetrical about the axis x-x. The 
axial stress a is constant through the thickness; it must 
be determined by conditions of axial equilibrium. Stress, 
deformation, and deflection equations can be written as 
follows : 




(*M a - CE « if teeif 



( 1 ) 




( 2 ) 
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Equations (4) and (5) give the outward displacement of 
inner surface and outer surface of thick-walled cylindrical 
vessel respectively. 

X. LONG THICK-WALLED CYLINDRICAL PRESSURE VESSELS 

UNDER INTERNAL PRESSURE 

Stress, deformation and deflection equations can be 
obtained from eqs. (1), (2), (3)j (4) and (5) in Section IX 
by simply putting p Q = 0, so as to obtain the following 
results. 

^ ft- _ feffe.'i R ’ 4 j_ 
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XI. LONG THICK-WALLED CYLINDRICAL PRESSURE VESSELS 

UNDER EXTERNAL PRESSURE 

Substituting = 0 in equations (1), (2), (3) s (^)> and 
(5) of Section IX we are led to the equations 

<r te'fa H ( j. 

02..ME, rtf + tr)2 (1) 

<s t _ _ fo.? % _ (p«y _ j_ 

lU-lK ;) 1 re„f-(E.;f lr, “" (2> 

(51') = - J . *. 

(P<o) 3 - ( ) a ° (at lnside surface) (3) 
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XII. THICK-WALLED SPHERICAL PRESSURE VESSELS 
UNDER INTERNAL AND EXTERNAL PRESSURE 

If the thickness of a spherical vessel is relatively 
large, the stress variation from inner to outer surface can 
not be neglected. Therefore, membrane theory can not be 
applied. For this case the spherical vessel can be analyzed 
by applying equilibrium equations of statics to an element 
also considering relations of geometric compatibility [*1,36]. 
Stresses, strains, etc., are functions only of radius and 
not of latitude or longitude on the sphere. 

Equations for stresses and deflections can be written as 
follows : 
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XIII. THICK-WALLED SPHERICAL PRESSURE VESSELS 



UNDER INTERNAL PRESSURE 



Stresses, deformation and deflection equations can be 
obtained from equations (1), (2), (3), (*0, and (5) in 
Section XII by putting p Q = 0 as follows: 
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(1) 



2 lrp[(B 10 ) 3 -(E s ,-f] 



( 2 ) 



K) 



max " 2 



(*».)* -lEs^ 



H- 



( 3 ) 



M. 






tv 



(it) 









3 V? 






k- 



( 5 ) 



28 



( 6 ) 



So 



1 j. , h a 




J = 2 P 





H- 



XIV. THICK-WALLED SPHERICAL PRESSURE VESSELS 
UNDER EXTERNAL PRESSURE 

t 

Substituting = 0 in equations (1), (2), (3), (4), 
and (5) of Section XII gives 



=• - 






A 






h 








1 I 

Jo 

4^ 

L ^1 


( r ) 3 1 






a 



♦ (P rf ) 

: — 7 1 °, 



(l) 



( 2 ) 



(3) 



29 



I '')' 4 [ tfs . f - IM 1 ] 



n= - - 

It 



(^) 



^s> = ^j_ * - j^' \ 3(i-v) • 



(?»)' 



r_E 



LV-tM' j 



Cl 



4" 



(5) 



.3 ^ \3 



^ li ) -. jj ? 4 . \ jv-ii-v). 

r = ? » 2E I iM*-«E, { y 



( 6 ) 



XV. ANALYSIS OF THIN-WALLED CYLINDRICAL SHELLS 
WITHOUT CLOSURE UNDER END LOADING 

This case is important in yielding coefficients 
(influence coefficients) required for dislocation analysis 
due to the effect of end closures or junctions with other 
shells, etc. The’ analysis of this case is given in Appendix 
1 (page 115 ). Only certain results are presented here in 
the body of the thesis. 

A. SEMI-INFINITE CASE 

The end deflection <5 ^ and rotation 0 ^ can be written 

c c 

in terms of end shear H-^ and end moment (see Fig. 1) 
as follows 
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pj_g # i _ End region of semi-infinite cylinder 
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B. FINITE CASE 

IP 1 

The deflections (6 C ,& c ) and rotations (0 C , © c ' 

ends 1 and 2 may be written in terms of the end 

(H 1 ,H 2 ) and moments (M 1 ,M 2 ) as follows (see Fig. 2) 
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) at 
shears 




Fig. 2 - Cylinder of finite length 
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In these equations the following abbreviations are used 

= Sintv fict 

Q = CO^Cv (l c t 

$ a = SIM. (3 c C 

^ - cos (3c f. 



The superscript 2 is not an exponent but, instead, indicates 
evaluation at the right end where x = A. 

XVI. THIN-WALLED CYLINDRICAL SHELLS WITH PLATE-HEAD 
CLOSURE UNDER INTERNAL AND EXTERNAL PRESSURE 

The analysis reported in this section and subsequent 
sections XVII through XXXI is for semi- infinite thin-wall 
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cylindrical shells which have various forms of closure at 
the accessible end. For a shell of finite length &, if the 
product 




is large (i.e., greater than 5.0, say) then the dislocation 
at one end has negligible effect on the conditions at the 
other end, and the analysis of a semi-infinite shell will 
provide satisfactory results. In shells which are not 
sufficiently long to justify this approximation, a more 
elaborate analysis is required. Brock [2] has treated the 
problem of cylinders composed of endwise juxtaposed short 
cylinders, each of constant thickness and properties, using 
matrix methods. The only type of closures explicitly consid- 
ered in [2] are (a) hemispherical and (b) ideally rigid; 
however the procedure can easily accommodate to other types 
of closure. 

Analysis of thin-walled semi-infinite cylindrical shells 
with flat head closure [7 , 15 , 23-a, 28 , 37 , 42 , 43a] , can be done 
by finding the deformations of each component separately, 
due to the pressure and arbitrary edge loads. The resultant 
deformations can be obtained by superposing the deformations 
due to the pressure and edge loads. Then the condition of 
continuity may be applied to provide enough equations to 
determine the edge loads necessary to produce continuous 
deformations (Fig. 1). 
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Fig. 1 Semi-infinite cylinder with flat plate closure 



At the junction (T) edge deflection and rotation 
equations for the plate can be written as follows; 
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The edge deflections and rotations of the cylindrical 
shell can be written in terms of loadings and the edge 
influence coefficients as follows 




r iH. 

“ £ c H, 



0 iM, 
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The edge influence coefficients are: 
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Conditions of continuity give two simultaneous equations 
which permit one to determine the edge loads as follows; 
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H ' f -gi — M i < — 5 — ( H--K) . 
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Once the edge loads are known, the stresses can be 
obtained easily by adding the contributions from the pressure 
and edge bending moment. 

XVII. THIN-WALLED CYLINDRICAL SHELLS WITH PLATE-HEAD 
CLOSURE UNDER INTERNAL PRESSURE 

Substituting p Q = 0 into the equations (1), (2), (3), 

(*0, (6), and (7) respectively of Section XVI gives; 
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( 3 ) 
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XVIII, THIN-WALLED CYLINDRICAL SHELLS WITH FLAT-HEAD 

CLOSURE UNDER EXTERNAL PRESSURE 

Substituting p. = 0 into the equations (1), (2), (3), 
(*0, (6) and (7) in Section XVI respectively 
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XIX. THIN-WALLED CYLINDRICAL SHELL WITH CONICAL-HEAD 
CLOSURE UNDER INTERNAL AND EXTERNAL PRESSURE 



The edge loads H^ and can be determined from compati- 
bility of displacement and rotation at the junction (Fig. 1). 
When the edge loads are known, the shear, circumferential and 

axial stresses can be determined at the junction [37, 42, 42a, 

. , <Pi-Po )d c 

43] . Each stress is divided by ^ (the circumferential 
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stress in a thin unrestrained cylinder under uniform 
pressure) to form a stress index, denoted by I. The stress 
index has two advantages. It reduces the magnitude of the 
results and enables one to tell immediately whether a stress 
is greater or less than the hoop stress in a thin cylinder. 
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Fig. 1 Semi-infinite cylinder with conical closure 

Some discrepency is to be expected between the theoret- 
ical results and the experimental results because of the 
effect of a transition zone a,b,c, (Fig. 1) where conical 
shell thickness may be considered as variable. 



The edge deflection and rotation equations for the 
cylindrical shell can be written in terms of edge influence 
coefficients as follows: 
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The edge influence coefficients are: 




(note that m = Vl2(l - v^) ) 

Similarly the edge deflection and rotation equations for 
the conical shell can be written in terms of edge influence 
as follows: 



41 






pM oi Mi ,p 

J/f ,, P 1 

A/1, + H, 



(fr-fcMf tfr-feJWcf in--fc)4 



(^0 



ZliJflL = _ 



O^ 1 0 iHi ,?i 

— *,i -^-«, + €U. 



( 5 ) 



The edge influence coefficients are: 
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(berj ^ + ( be i x 

C = $o be $'jj 0 + beij b e ij[ ) 

G * ( bir S *jo f f ljL$of 

The formulas for the Bessel (type) functions ber and bei 
can be found from [24-a]. 

It is not necessary to calculate the edge influence 
coefficients for the conical shell; numerical values have 
been tabulated in [43], including proper algebraic signs so 
that the following formulas are correct 

Conditions of continuity give the following equations 
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Simultaneous solution of equations (7) and (8) gives the 
edge loads. 
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Stress index equations for the cylindrical shell can be 



written as follows; 
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Similarly stress index equations for the conical shell 
can be written as follows: 
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These stress indices are tabulated in [^3]. 



XX. THIN-WALLED CYLINDRICAL SHELLS WITH CONICAL 
CLOSURE UNDER INTERNAL PRESSURE 

Substituting p Q = 0 into the equations (7), (8), (9), 
(10) j (11), (12), (13) and (lM) of Section XIX gives 
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XXI. THIN-WALLED CYLINDRICAL SHELLS WITH CONICAL-HEAD 



CLOSURE UNDER EXTERNAL PRESSURE 



Substituting p.^ = 0 into the equations (7), (8), (9) s 
(10), (11), (12), (13), and (1*1) of Section XIX gives 
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XXII. TWO ADJACENT CYLINDRICAL OR HEMISPHERICAL SHELLS 
OR CYLINDRICAL SHELL WITH HEMISPHERICAL HEAD 
CLOSURE UNDER INTERNAL AND EXTERNAL PRESSURE 

The case of two abutting cylinders is contained in what 
follows by taking the parameters n^ = n 2 = 2 . The case of 
abutting hemispherical shells is given by n^ = n 2 - 1. The 
case of a cylindrical shell closed by a hemispherical shell 
is given by taking n^ = 2n 2 = 2 or n 2 = 2n^ = 2. That all 
these cases can be treated by the same method, i.e., by 
assuming, in effect, that the sphere behaves near its edge 
as if it were a cylinder of equal radius, results from the 
fact that the edge dislocation effects are confined to a 
narrow band near the edge, and in this band the spherical 
shell is actually very nearly a cylinder. The difference 
between the several cases thus arises solely in the differ- 
ence between the change in radius due to pressurization in 
a cylinder and a sphere [4,15,20,25,33,37,42].- 
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CytmcJer (*1, = ^) 
or Sphere i 



Ci^iiAdsc C fla.-3.) 

or Sphere ( tOj. = i ) 



I 

Fig. 1 Abutting cylindrical or hemispherical shell 

elements 



We have 








( 2 ) 
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and thus obtain 




which can be solved to yield M and H. 

If the thickness and the material of each component 
is the same, then obviously 



M ^ 
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XXIII. THIN-WALLED CYLINDRICAL SHELLS WITH HEMISPHERICAL- 



HEAD CLOSURE, WITH SAME MATERIAL AND SAME THICK- 
NESS UNDER INTERNAL PRESSURE 



The formulas of Section XXII give 



M s o 




(i) 



( 2 ) 



XXIV. THIN-WALLED CYLINDRICAL SHELLS WITH HEMISPHERICAL- 
HEAD CLOSURE, WITH SAME MATERIAL AND SAME THICK- 
NESS UNDER EXTERNAL PRESSURE 



The formulas of Section XXII give 



M - o 




(i) 



( 2 ) 



XXV. THIN-WALLED CYLINDRICAL SHELLS WITH SPHERICAL 
HEAD CLOSURE UNDER INTERNAL AND EXTERNAL PRESSURE 

Since the loading is axisymmetrical , the deflection and 
rotation equations of both cylindrical and spherical shells 
at the junction can be written in terms of edge influence 
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coefficients [13 , 18 , 37 , ^2 ] . Then the conditions of continu- 
ity for the junction (Fig. 1) give simultaneous equations 
for edge shear force and bending moment. If the angle <f> = 90° 
so that the spherically curved closure becomes a hemisphere, 
the formulas herein simplify greatly and correspond to those 
already given in Section XXII. 




Fig. 1 Semi-infinite cylinder with spherically shaped 

end closure 



The edge deflection and rotation for the cylindrical 
shell can be written in terms of edge influence coefficients 
as follows 
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The edge influence coefficients are: 
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The edge deflection and rotation for the spherical shell 
can be written- in terms of edge influence coefficients as 
follows : 
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The edge influence coefficients are: 
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Conditions of continuity give the following simultaneous 
equations 
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XXVI. THIN-WALLED CYLINDRICAL SHELLS WITH SPHERICAL-HEAD 



CLOSURE UNDER INTERNAL PRESSURE 

Substituting p o = 0 into the equations (1), ( 2 ), (*0, 
(5)> (7), and (8) in Section XXV give: 
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XXVII . THIN-WALLED CYLINDRICAL SHELLS WITH SPHERICAL-HEAD 

CLOSURE UNDER EXTERNAL PRESSURE 



Substituting = 0 into the equations (1), (2), (4), 
(5), (7), and (8) in Section XXV give*. 
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XXVIII. THIN-WALLED CYLINDRICAL SHELLS WITH TORI SPHERIC AL- 
HEAD CLOSURE UNDER INTERNAL AND EXTERNAL PRESSURE 



The analysis of thin-walled cylindrical shell with 
torispherical head closure (Fig. 1) can be accomplished by 
knowing the behaviors of three separate parts. Edge 
deflection and rotation equations can be written in terms 
of edge influence coefficients and then the condition of 
continuity may be applied [11,12,13,15 J 37>^2]. 





Fig. 1 Semi-infinite cylinder with closure that is partly 
toroidal and partly spherical (i.e., 
"torispherical. " ) 



For the cylindrical shell at junction 1 , the deflec- 
tion and rotation can be written in terms of edge influence 
coefficients as follows: 
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The edge influence coefficients may be taken from 
equation (3) in Section XXV or alternatively can be written 
as follows: 
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For the toroidal shell at junction (T) the deflection 
and rotation can be written in terms of edge influence 
coefficients as follows: 

(4) 

(5) 
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Similarly for junction ( 2 ): 




( 6 ) 



( 7 ) 



In the toroidal shell it is not permissible to assume 
that its two edges are so remote from one another that their 
mutual interaction may be neglected; accordingly the formulas 
for the edge influence coefficients are rather complicated 
and will not be reported here. Necessary information may be 
found in [11,12], However, as before, the cylinder may be 
assumed to be of semi-infinite length and the spherical 
portion is assumed to subtend a substantial angle, (i.e.. 



For the spherical shell at junction , the deflection 
and the rotation can be written in terms of the edge 
influence coefficients as follows: 
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The edge influence coefficients can be written as follows 
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Conditions of continuity give four simultaneous equations 
to solve, as follows: 
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( 12 ) 
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Simultaneous solution of equations (11), (12), (13) and 
(1*0 gives the edge loads. 



XXIX. THIN-WALLED CYLINDRICAL SHELLS WITH 
TORISPHERICAL-HEAD CLOSURE UNDER 
INTERNAL PRESSURE 



Substituting p Q = 0 into the appropriate equations in 
Section XXVIII gives 
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XXX. THIN-WALLED CYLINDRICAL SHELLS WITH TORISPHERICAL 
HEAD CLOSURE UNDER EXTERNAL PRESSURE 

Substituting = 0 into the appropriate equation in 
Section XXVIII gives 
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( 5 ) 
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XXXI . THIN-WALLED CYLINDRICAL SHELLS WITH 
ELLIPSOID-TORUS TRANSITION 
UNDER INTERNAL PRESSURE 



Analysis can be performed by writing the deflection and 
rotation equations for each junction in terms of the edge 
influence coefficients. Then the condition of continuity 
can be applied to obtain a set of simultaneous equations 
[13]. Two cylindrical shells connected by an ellipsoidal- 
torus transition have been taken as a typical example 
(Fig. 1). 
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Fig. 1 Ellipsoid-torus transition between 

coaxial cylinders 
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For upper cylindrical shell the deflection and rotation 
equations at junction (T) , can be written in terms of edge 
influence coefficients as follows: 




= S 



Mi 






r^/H/ 

<>C H, 






-t 




( 1 ) 



(2) 



For the torus the deflections and rotations at junctions 
© and © respectively, can be written in terms of edge 
influence coefficients as follows: 
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For ellipsoidal shell the deflection and rotation at 
junction ( 2 ) and (3) respectively, can be written in terms 
of edge influence coefficients as follows: 
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( 7 ) 



( 8 ) 



( 9 ) 



(10) 



For the lower cylindrical shell the deflection and 
rotation at junction (3) can be written in terms of the edge 
influence coefficients as follows: 
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Conditions of continuity give the following simultaneous 
equations : 

pi/Ui, r 1 ^ /nlHi-iH/. / fs/Pf 0 ip< \ u _ (13) 



/0 aM2 n 3Ma. n 2^3 -jj. 

it -5 e J*-Se /«,♦£ VfA 34 )* n 3 * ,R jP 

3 *«,*<* ie 



Je (15) 



Vv«Vj* ,- e;Ve;V • » 






Af? t( ®T- 9^ 3 )M 3 + 9g \ t l Q^- of 2 ) H 3 + Pj^ ^ (l8) 
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Tabulated numerical values of these influence coefficients 
may be found in [13] and in references cited in [13]. Close 
attention should be given to the sign conventions employed 



Simultaneous solution of equations (13), (1*0, (15), (16), 
(17) and (18) gives the edge loads. 

If the toroidal and ellipsoidal shell segments are "long", 
that is if their meridional arc lengths are several times 
the local value of an attenuation length, then the "cross- 
influence" coefficients (reflecting the effect of conditions 
at one edge upon those at the other) are quite small compared 
to the "self-influence" coefficients and may be neglected. 

This affords a considerable simplification, as follows: 



in all these references, including [12]. 




( 22 ) 



( 23 ) 



( 21 ) 



(19) 



( 20 ) 
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( 24 ) 



(25) 



( 26 ) 



Then the simultaneous equations (13), (14), (15), (16), 
(17), and (18) become: 




(27) 




( 28 ) 





(29) 



( 30 ) 
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(31) 



tK h -g* h ) fe. = 



e J r z - o 



(<£ M - ef> 3 1 ( of 3 - <ef s j + ©e\> (32) 



These equations are not only simplified in appearance 
but are now uncoupled pairwise so that (27) and (30) yield 
M 1 and H-^ (28) and (31) yield M 2 and H 2> and (29) and (32) 
yield and H^. 

The reader may observe that in this section an ellip- 
soidal segment is involved although we have not elsewhere 
given a treatment for a full ellipsoidal closure. Data for 
influence coefficients of so-called "open crown" ellipsoidal 
shells is more readily available than for complete semi- 
ellipsoidal closures. See Galletly [13] and references 
cited therein. Kraus et al. [22] have used a digital computer 
to analyze thin-walled cylindrical vessels with full semi- 
elliptical heads but the results are not in the form of 
influence coefficients. Instead, coefficients of (stress 5 
pressure) are tabulated for useful ranges of several geo- 
metric parameters describing the configuration. It is not 
in the spirit of this thesis to reproduce this tabular data 
and the reader is referred to [22]. 
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XXXII. AXIAL STRESS IN THIN-WALLED CYLINDRICAL SHELLS 



UNDER INTERNAL AND EXTERNAL PRESSURE 
TOGETHER WITH AXIALLY APPLIED FORCE 

Uniform distribution of the axial stress on the cylinder 
cross section has been assumed. The axial stress is then a 
function of the geometry of the cylinder [ 6 , 32 , 33 , 3 ^]* 

Axial stress can be written as follows: 




( 1 ) 



(Note that F = 2 ttR q is externally applied axial tensile 
force) 



If only internal pressure is applied equation (1) becomes 




( 2 ) 



If only external pressure is applied equation (1) becomes 




* A 



3 - 



( 3 ) 
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XXXIII. AXIAL STRESS IN CLOSED THICK-WALLED CYLINDRICAL 



PRESSURE VESSELS UNDER INTERNAL AND 
EXTERNAL PRESSURE 

Axial stress is uniformaly distributed on the cylinder 
cross section [6,32,33»3^] except near the junction of 
cylinder and closure. If there is no axial loading other 
than that due to pressure the axial stress can be written as 
follows : 



(1> 

If only internal pressure is applied equation (1) becomes 




( 2 ) 



If only external pressure is applied equation (1) becomes 



<n„f % 
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XXXIV. 



AXIAL STRESS IN THICK-WALLED CYLINDRICAL PRESSURE 



VESSELS UNDER INTERNAL AND EXTERNAL PRESSURE 
TOGETHER WITH AXIALLY APPLIED FORCE 



In this case the axial stress equations can be written 
as follows : 



^ - ( ^ C£ >) 'fe + IT 

cm*- c Zdf 



CD 



where F is the externally applied axial tensile force. 

If only internal pressure is applied equation (1) 
becomes 







( 2 ) 



If only external pressure is applied equation (1) 
becomes 







( 3 ) 



XXXV. BUCKLING OP THIN-WALLED CYLINDRICAL SHELLS WITH 



OPEN ENDS UNDER EXTERNAL PRESSURE 



The analysis is discussed in [4,8 ,15>33>35>37>42] . 

Only the results will be given here. It is assumed that 
there is no restraint against axial deformation and that 
the only restraint at the ends is agains radial deformation. 
The general solution can be written as follows 



where n is an integer which must be selected (by trial or 
otherwise) so as to minimize the value of this expression. 
The value of n, however, must be not less than 2. 

If the ratio d /% is very small, the expression simpli- 



fies considerably and the value n = 2 is then seen to be 
appropriate so that. 



The preceding equation presumes elastic action. The 
following formula Is for elastic action with no strain 
hardening [33]. 




( 1 ) 




4 0-vfKe^ 



( 2 ) 




( 3 ) 
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XXXVI. BUCKLING OF THIN-WALLED CYLINDRICAL SHELLS WITH 



CLOSED ENDS UNDER EXTERNAL PRESSURE 
General buckling pressure can be written as follows [4, 

333: 




n is an integer which must be selected (by trial or other- 
wise) so as to minimize the value of this expression. 

The end constraints are as in Section XXXV except that 

2 

an axially compressive force P = p(d ) /4 is also present. 

c o 

Formula (1) is cumbersome to use, so that the following 
approximate formula in which n, the number of lobes of the 
buckled shape does not appear, has also been employed in 
this case 



6 
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XXXVII. SYMMETRICAL BUCKLING OP THIN-WALLED CYLINDRICAL 



SHELLS UNDER ACTION OF UNIFORM AXIAL COMPRESSION 



The configuration of the system is shown in Fig. 1. 




Fig. 1 Buckling of thin-walled cylinder due to 
due to uniform axial compression 



The critical axial compressive stress is 






xc 



m 



(r? trf 2 etc 



(pir^Ecte. 13. ti-WJff) ! 



(i) 



where n is an integer (n _> 1) which must be chosen so as to 
minimize this expression; n represents the number of half 
sine waves into which the cylinder buckles. The half wave- 
length is &/n. The value of n is given approximately by 

Wc. 

TT 



n = 



( 2 ) 
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and if this value is substituted in Equation (1) the 
following approximation is obtained 




■Ec. 



( 3 ) 



\] 3 



The case n = 3 is shown in Fig. 1. 

When the shell is not very thin and buckling occurs at a 
stress that is beyond the proportional limit the critical 
load cart be obtained by introducing the tangent modulus, E^, 
instead of the elastic modulus E at an appropriate place in 
the analysis. Then the approximate result becomes 



If the thin-walled cylindrical shell is short, i.e., if 



the critical buckling stress can be written as follows: 





(4) 







IT 




g- _ inf B (tcf 



(5) 
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XXXVIII. BUCKLING OF LONG THIN-WALLED CYLINDRICAL 



' SHELLS SUBJECTED TO TORSION 

Thin cylindrical shells may buckle as a whole without 
flattening of the circular cross section, or may buckle 
locally with modification of cross section. In the latter 
case if the shell is long enough, the critical moment does 
not depend essentially upon the conditions of end' restraint , 
and is given by [35] 



M 



cr - 




(l) 



From the simple relation 



Mcr — 2 n (2 C ) "Ccr 

we obtain 




( 2 ) 



( 3 ) 



For shorter cylinders, the modification of cross section 
involves several (i.e., more than two) lobes and the 
formulas in Section XXXIX and XL should be used. . 
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However, the cylinder may also buckle as a whole. If 
no lateral restraining moments are applied at the ends, the 
corresponding critical moment is given by 




a(frf C2 c f-h c E 
£ 



(*0 



which may be smaller than the result given by Eq. (1) of 
the present section or the analysis of Sections XXXIX or XL, 
whichever may otherwise be applicable. Similar results are 
not available for other than laterally momentless end 
conditions . 



XXXIX. BUCKLING OF THIN-WALLED SHORT CYLINDRICAL 

SHELLS WITH CLAMPED EDGES SUBJECTED TO TORSION 

If [35], 



t < n.a 




(1) 



then the critical shear stress can be obtained from 



(-»£) 

E 




1 
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3M4 



( 2 ) 



and the corresponding value of may be obtained from 
Equation 2 of Section XXXVII. Cf. also Equation k of 
Section XXXVIII. 
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XL. BUCKLING OF THIN-WALLED SHORT CYLINDRICAL SHELLS 



WITH SIMPLY SUPPORTED EDGES SUBJECTED TO TORSION 

If C35], 




( 1 ) 



then the critical shear stress can be obtained from 







1 + 0 . 0257 ; 







( 2 ) 



and the corresponding value of M may be obtained from 

O X 

Equation 2 of Section XXXVIII. Cf. also Equation 4 of 
Section XXXVIII. 



XLI . BUCKLING OF UNIFORMLY COMPRESSED THIN-WALLED 
SPHERICAL SHELLS UNDER EXTERNAL PRESSURE 



The theoretical critical buckling pressure and stress 



are [35] 









(i) 




£ 7l s 



( 2 ) 



81 



However, externally pressurized thin spherical shells 
seem to be quite sensitive to the presence of slight 
imperfections and to small disturbances during loading, and 
observed values of critical pressure are usually much smaller 
than given by Eq. (1). Accordingly, the following empirical 
formula is offered [35] 




( 3 ) 



which agrees with experiments for 



400 < 



- t 



< 2000 and 20° < 4>° < 60 ° 



c 



See Fig. 1. 




Fig. 1 Spherical dome under external pressure 



82 



XLII. THIN-WALLED CYLINDRICAL SHELLS UNDER CENTRIFUGAL 



FORCE FIELD 

When the cylindrical shell is rotated about its axis 
with constant angular velocity w (radians per second) uniform 
radial loading is produced [33»3^I]. By considering radial 
equilibrium of an element cut from a cylinder and assuming 
that there is no axial load or stress, one obtains 

^ ( 1 ) 



The radial deflection is 

£ = -£ (wf ( 2 ) 

and the axial unit strain is 

(UjJ (3) 

If a x 7 £ 0, the following formulas should be used in 
place of (2) and (3) 
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XLIII. THICK-V/ALLED CYLINDRICAL PRESSURE VESSELS 



UNDER CENTRIFUGAL FORCE FIELD 



This case is similar to but not the same as the well 
known case of a rotating thin disk. The latter is a case 
of plane stress whereas we are here concerned with a case 
of plane strain. 

Specializing first to the case e = 0, it may be shown 

A 

that the formulas for a and a. in the case of a thin disk 

2 

may be used if the angular velocity squared, to , appearing 
in the formulas, is replaced by a modified value, namely 

— 2 of! - 2 v ) 

(to) = (to) 4- . The resulting formulas are 

(1- v) 2 






( 1 ) 



< 4 = 3 
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3* V 



where 



^ _ (3^)0 +^(1-2^) P(^) a 
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(3) 



( 4 ) 



The formula for a comes from the condition e - 0; 

X A 

namely, a x = v(a r + cr t ). 



\2 



(Sr) m3x = B( at r = \J2 U V. 



Cl 



( 5 ) 
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at r-B Cl ' (6) 



The radial deformation is 





Thus, also 




2(hV) 

(i+w)(i-av// E 
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( 8 ) 
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( 9 ) 
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If e x / 0, we must superpose the stress system 




(10 a,b,c) 



and the corresponding additional deformation 



It = - r £ x 



(ii) 



Thus , we have 




to be added to the results given in Eqs. (8) and (9) 
respectively . 

XLIV. THIN-WALLED CYLINDRICAL SHELLS WITH RADIAL 

TEMPERATURE GRADIENT 

The assumption of thin wall leads to the approximation 




0 



( 1 ) 



It is convenient to define an average temperature 

Ze, 



T = 





( 2 ) 
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Then 




( 3 ) 




V' I ~ \lr) +E^ 



( 4 ) 



This corresponds to an axial tensile force. 



T - tA ( £ n - oil ) 



( 5 ) 



where A is the annular cross section area. If there is no 
axial force e = aT. 

XLV. THIN-WALLED LONG CYLINDRICAL SHELLS WITH 



It is assumed that the temperature at any point is a 
function of axial position only. There is no variation 
through the thickness of the wall. In particular, also, the 
following results apply to the case where a segment of thin 
cylinder with constant linear temperature gradient is 
joined to a segment of the same pipe having constant temper- 
ature. A more general situation can be treated by integrating 
a fourth order differential equation [3*1]. 

The general equation for symmetrical deformation of long 
cylindrical shell can be used for analysis [5 ,21 , 33 , 3*1 ] • 



TEMPERATURE GRADIENT IN AXIAL DIRECTION 
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Denoting by T the Increase of the temperature of the shell 
from a certain uniform initial temperature and assuming that 
the shell is divided into infinitely thin rings by planes 
perpendicular to the x axis, the radial expansion of the 
rings due to the temperature changes is aR T. This expansion 
can be eliminated by applying external pressure of an 
intensity which brings the shell to its initial 

diameter. Assume that the length b (Fig. 1) of the thin- 
walled cylindrical shell is long and the part of the cylin- 
drical shell to the right of the cross section 1-1 has 
constant temperature T^, and to the left has a temperature 
which decreases linearly to the value (T + T^) at the end 
x = b . 



t 




Fig. 1 Cylinder having linearly varying temperature 

adjacent to section having constant temperature 
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Then the temperature increase, as a function of x, is 



T « 



To x 
b 



(i) 



Using this in the general differential equation and 
applying boundary conditions, the equations for bending 
moment and maximum stress at section 1-1 can be written, 
respectively, as follows 



M b = ft c b c 'Ll* To 

2b 



( 2 ) 




max 



Oh- = o.353 
0 tcf b 



( 3 ) 



It has been assumed that the length b >> R t . Other- 
wise a correction to the moment equation becomes necessary 
and can be written as follows; 



-2(3cb[ 

4 a WM b e 



^COS(PctJrSh(0cW)L < 3 / v| 1>e 



This moment correction should be applied to the moment 
given by equation (2) to get the correct moment,- after which 



6 (M b + AfA h ) 






a 



(5) 



XLVI . THIN-WALLED CYLINDRICAL SHELLS UNDER INTERNAL 



PRESSURE, AXIAL TENSILE FORCE AND NONUNIFORM 
TEMPERATURE CHANGE 



This section is adapted from the treatment in [2,24]. 

By considering the equilibrium of the element shown in Pig. 1, 





x 



Fig. 1 Element of thin cylindrical shell showing 

dimensions and notations 



the following equations can be written 

% * ^ 6 * cLr (1 > 

±, 

~ a 
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■t C 







Also the moments can be written in terms of 
as follows : 




By considering the geometry shown in Fig. 2 
ing equations can be written 
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( 3 ) 



the stresses 



(4 a ,b ) 



the follow- 



( 5 ) 



( 6 ) 
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Fig. 2 Drawing showing geometry of deformed shell 

where subscripts m and b denote the membrane component and 
bending component respectively. It has been assumed that 
the strain variation is linear through the wall. Also 
T = T(x) denotes the local average temperature through the 
wall and AT = AT(x) denotes the local excess of outside 
wall temperature, the variation being assumed linear. 

We also have 

- oi AT 

k = — - th *~ a 

? ' ±c 

a 

where k is longitudinal curvature. 
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Making the customary approximation 



_ OS* 

dX 1 'to 



Since the cross section preserves its circularity we 
can write 




( 8 ) 



From equations '( 6 ) and (8), by comparing like powers of 
r, it may be seen that 




U AT 
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( 9 ) 
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( 10 ) 



By considering the stress-strain relations, the following 
equations can be written 




L 



(ID 
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( 12 ) 





(£ m t + W£ m> )-( J _ r J(j t ttV^£ M ) 



The radial stress has been neglected in eqs . (11) and 

( 12 ). 

Manipulation of these equations leads to the differential 
equation 




For a shell with constant properties along the length, 
this reduces to 

-&t ~ 4- ■ — + T (lio 

This fourth order differential equation must be solved 
subject to appropriate end conditions which will be discussed 
later in this section. Once u is known as a function of x, 
stresses may be found from the formulas 
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(16) 




where 




d *IA (HVO d .AT 

*t c 



(17) 



The end conditions with which we may be concerned relate 
either to a condition at an accessible edge or to a condition 
at an infinite distance. At an accessible edge, any two of 
the following four conditions must be specified: 

a) end displacement u at x = 0 

b) end slope ^ at x = 0 

c) end moment (see Eq. 17) at x = 0 

. dM x 

d) end shear at x = 0. 

At an infinite distance, we must have vanishing shear 
and vanishing slope, with radial displacement satisfying 



and moment satisfying 



Dc ( 1+ K) 




(19) 



A total of four such conditions is required to specify 
a unique solution to equation (14), two from each end, where 
the word' end means either an accessible edge or an infinitely 
distant section. 

Since such a wide variety of conditions is represented 
by this analysis, no attempt will be made to give specific 
applications, except that we may note that if the right hand 
side of eq. (l^J) is constant, we have 



where 3 and Q are constants. The general solution of this 
equation is 



where C^, C2, and are arbitrary constants to be deter- 
mined by the end conditions and h^(3x) are functions defined 
and discussed in Appendix 2. 




( 20 ) 
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XLVII . THICK-WALLED CYLINDRICAL PRESSURE VESSELS WITH 



' RADIAL TEMPERATURE GRADIENT 



Thermal stress analysis of thick-walled cylindrical shells 
is extensively reported in pressure vessel literature [3,5,6, 
19 , 21 , 33 , 34 , 36 , 42 ], 

Analysis can be accomplished by taking an element and 
then applying to it equilibrium equations of statics, stress- 
strain relations, and compatibility. It has been assumed 
that the temperature is symmetrical with respect to the axis 
of the cylinder and constant along this axis. Therefore the 
deformation of the cylinder is symmetrical about its axis. 
Plane strain has been considered. 

Equations for radial stress, circumferential stress, and 
deformation can be written, respectively, as follows: 
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If the axial force F is zero, then both ends of the 

X 

cylinder are free of traction. 

Equations for 0 , a . , a and u can be written, respec- 
tively, as follows: 
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(7) 



Note that equations (4) and (5) are the same as with 

equations (1) and (2) respectively, because 0 and 0 . are 

r* u 

not dependent on 0 . 
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If e =0, then the ends of the cylinder are not free of 

X 

traction. Equations for a , a. , a and u can be written, 

x T/ X 

respectively, as follows: 
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Again the stresses a v and are not changed. 

If the distribution of temperature T(r) as a function 
of radial position is known, the integrals in all above 
equations can be calculated; then the stresses and deformation 
may be obtained at any desired position. 
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XLVIII. THICK-WALLED CYLINDRICAL PRESSURE VESSELS 
WITH LOGARITHMIC RADIAL TEMPERATURE GRADIENT 

Steady state radial heat flow through the walls of a 
pressure vessel, an obviously important special case, leads 
to the so-called logarithmic gradient given by 




where T^ represents the temperature of the inside surface and 
where, without essential loss of generality, the temperature 
of the outside surface is taken to be zero. 

Using this particular temperature distribution in the 
formulas of the preceeding section leads to those given in 
the present section. 

If the total axial force is zero and the ends of the 
cylinder are free of traction, equations for a , o. , a and 
u can be written, respectively, as f ollows : 




100 






(4) 






The maximum circumferential stress, which is of greatest 
interest, occurs at the inner surface [1-a]. 




( 6 ) 



If e x is zero and the ends of cylinder are not free of 
traction, equations for o , a, , a and u can be written, 

x o X 

respectively, as follows: 



- The same as 


(2) 


(7) 


= The same as 


(3) 


(8) 
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( 9 ) 
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Maximum circumferential stress, which occurs at the 
inner surface, is given by Equation (6) above. 



( 10 ) 
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XLVIX. THICK-WALLED SPHERICAL PRESSURE VESSELS WITH 



STEADY STATE TEMPERATURE GRADIENT 

The well known [36] radial and circumferential stresses and 
deformation equations can be written, respectively, as 
follows : 





r 4- 







(l) 



( 2 ) 



(3) 



In the case of steady state heat flow the temperature at 
any distance from the center can be written as 




where we take T(R ) = 0. 

c o 
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Substituting equation (4) into the equations (1), (2), 
and (3) gives: 




J. t Tj* 

0-v ) 





(5) 
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L. THICK-WALLED TWO-LAYER CYLINDRICAL PRESSURE VESSELS 
WITH DIFFERENT MATERIALS UNDER SHRINK-FIT TOGETHER 
WITH INTERNAL AND EXTERNAL PRESSURE 

Although the general case of multilayered vessels having 
n layers is given in Appendix 3> the analysis there involves 
some manipulation. Since the two-layered vessel represents 
the most important application, formulas for this particular 
case are presented here. (See also discussion at end of 
Appendix 3 . ) 
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Pig. 1 Dimensions and notation for two-layered vessel 



The interface pressure can be written as follows 








(i) 
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For the inner cylinder, equations for stresses and 



deformation can be written as follows: 




For the outer, cylinder equations for stresses and 
deformation can be written as follows: 
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LI. THICK-WALLED TWO-LAYER CYLINDRICAL PRESSURE VESSELS WITH 
SAME MATERIALS UNDER SHRINK-FIT TOGETHER WITH INTERNAL 

AND EXTERNAL PRESSURE 



The interface pressure can be written as follows: 
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( 1 ) 



For the inner cylinder equations for stresses and 
deformation can be written as follows: 
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( 4 ) 



For outer cylinder equations for stresses and deformation 
can be written as follows : 
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Equations for the cases with internal and external 
pressure only can be obtained by putting p Q = 0 and p^ = 0, 
respectively, into the equations given in Section L and in 
this section. 
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LII. OPTIMAL DESIGN OF MULTILAYER PRESSURE VESSELS 



OF SAME MATERIAL 



Before presenting the sequence of formulas, a remark is 
in order about what is meant by "optimizing" in this case. 

It means simply that the overall volume or weight of the 
metallic material is minimized, subject to the conditions 
that the inside radius is the given number R c ^, that the 
internal pressure is the given number p^, and that in the 
loaded condition each of the shells is stressed up to the 
allowable limit, but no higher (Fig. 1) 

It has been presumed that the design pressure p^, has 
been specified, and that the material selection has provided 
values of E, Young’s modulus, of v, Poisson's ratio, and of 
, the allowable tensile stress. It should also be 
explicitly stated that all shells are presumed to be composed 
of the same material, or, at least, that the values of S^., E, 
and v are the same for all shells. Also, one should be aware 
that the following formulas yield results which are precisely 
valid only at some distance from any "end effects". The case 
of plane strain is postulated [1]. 

Equations applicable to the choice of the maximum 
shearing stress theory of failure will be denoted as (MSS) 
and equations applicable to the choice of the maximum normal 
stress theory of failure will be denoted (MNS). 

If the MNS theory is used, S^ is simply the allowable 
normal or tensile stress; if (MSS) theory is used, S^ is 
twice the allowable shearing stress. 
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Note: Symbols a, S, < t 




Fig. 1 Dimensions and notation for multi-layer vessel. 



The sequence of design formulas can be written as follows: 




i- 

n. 

- (itf)t 



(i) 



(2) MSS 



(2) MNS 
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(b -=■ 



(3) MSS 



& = J— 

UM 

K = 



(3) MNS 



( 4 ) 



where K is the ratio of metallic volume to the internal 
volume, and thus is a measure of the "optimality" of the 
design. For both theories of failure, optimum design 
results in shell radii which are in geometric sequence. 



(P) T 



(5) 



The radial compression which exist in the fully assembled 
and pressurized condition between i^ h and (i+l)^* 1 shells is 




(6) MSS 



(6) MNS 
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The desired radial interferences are given by 




j ( i- K % 



*i 



Co 






(7) MSS 



£. = (7> MNS 

e 

These radial interferences are the amount by which the 

4 - 

outer radius of the i shell exceeds the inner radius of the 

t" Vi 

(i+1) shell when the shells have been manufactured and are 
at the same temperature and before the assembly process has 
started. 

The sequence of formulas for stress can be written as 
follows . 

These formulas relate to the problem of calculating 
stresses at various points of the assembled vessel in the 
unloaded and in the pressurized condition. 

Hoop stress: 





(8) MSS 



(8) MNS 
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(9) MSS 




U( (b-l ) 

Axial stresses: 



(9) MNS 
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(10) MNS 
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(11) MSS 







(11) MNS 



It should be noted that the formula (10) MNS gives only 

n i 

the single value a ; and to get other values of a for the 

X A 

(MNS) case, formula (11) MNS should be used. 

To complete the calculations of stresses, the formulas 
will be exhibited which show the changes in going from the 
unpressurized to the fully pressurized condition. 
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APPENDIX 1 



SEMI-INFINITE THIN-WALLED CYLINDRICAL 
SHELLS WITHOUT CLOSURE UNDER END LOADING 





Fig. 1 End region of semi-infinite cylinder 

The solution of the governing fourth order differential 
equation, which also satisfies the requirement that the 
radial displacement, u , and its derivatives vanish as x-* 00 , 
may be written as (cf. Section XL VI) 

(A — A I j^i (ficx) f A3 ^2 ( 1 ) 

where A-^ and A2 are arbitrary constants and f-^(3 c s) and 
^2^c X ^ are an ^ two the four functions 
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PcX -Pc A 

A(W e cps(^cX) JfPex)* c sm(^x) 

( 2 ) 

-^(M) - A((3 C *} +I((3 c x) (Pt*) *-A-(PcX)-X(/3 t xj 



which have the following values for zero argument 



_/L(o)= -Q<q) - YloJ =( 5/p) = 0 



( 3 ) 



and which have the schedule of derivatives shown in Table-1. 



OPERATION FUNCTION 
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Table-1 SCHEDULE OF DERIVATIVES FOR FUNCTIONS (2) 
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A convenient selection gives (u) in terms of the edge 
shear and moment M^. (See Fig. 1), viz. 



This gives the following evaluations of (u) and its 
derivatives at the accessible end x = 0 



, __ 

in- _ _ til-*- a$ c M { 

" 2(Uc f D 0 

dhi /Mi 

at x = 0 (5) 

__ ti. 



These results provide the influence coefficients frequent- 
ly encountered in the body herof; e.g. cf. Equation (3) of 
Section XXV. 




a i (3 c ) 3 v c 
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APPENDIX 2 



FINITE THIN-WALLED CYLINDRICAL SHELLS 
WITHOUT CLOSURE UNDER END LOADING 

In this appendix we discuss the behavior of finite thin- 
walled cylinder under edge (end) loadings and obtain the 
appropriate corresponding influence coefficients. 




Fig. 1 Finite cylindrical shell 

The functions 

= cask. (£cx) C3S(ficX) 

k 2 ((Z c X)= j[c3 (flcXj CsS (fax) j 

A 4 f(i c x) =A[ct.s-ei((i<5<)^w((3cx)-srnk(^) ceslhxj] 
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provide four linearly independent solutions to the homo- 
geneous form of Equation (2) of Section XLVI , the solution 
of which can thus be written 

^ — A , ^ t "t" (2) 

These functions have the following values for zero 
argument 

t 

A*®)- I \lo) *4i z (0) = 4 4 (9j ^ o ( 3 ) 



and their derivatives are shown in Table 1 
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-4 (&:)%(&*) 







Table-1 SCHEDULE OF DERIVATIVES OF THE FUNCTIONS h ± 
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Thus we obtain. 



dx = ' t(3 ^' ^ 



(4) 



d ? M 



d *» ^M(MJ t ft ») < s > 




From Equations (3) it is easy to obtain the evaluations 
at the left end (indicated by superscript 1) 



£,' = a, 

9l = ( 3 . &> 

A4, 

( (M* D c " 



(7) 

( 8 ) 

(9) 



-H i 



(10) 
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Similarly, equations (2), (4), (5), and (6) can be 
rewritten for the end 2 as follows: 







£ c = Q,it t -f ^2 * ^3 As + ^4^-4 

= -4MA* t0^ 3 ^t(3 c <?X 



i; -- - £ 

A T^c 



( 11 ) 



( 12 ) 



(13) 



(14) 



where 

kkitd) ^A^d) 

(Note use of superscript 2 to indicate evaluation at x = Z.) 

Equations (9) and (10) give coefficients a^ and a^ in 
terms of the end loadings. Thus, the only unknown quantities 
in equations (13) and (14) are the coefficients a^ and a2, 
and simultaneous solution of (13) and (14) gives 



t--W) 



(15) 
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(16) 





\-k%2 +40u] + W, ] + V -+42^2 



(17) 



We emphasize again the use of superscript 2 to indicate 
evaluation at the right end; exponents appear outside 
parentheses . 

These expressions for a^ and a.^ could be rewritten in 
terms of hyperbolic and trigonometric functions of (3 C &) 
and then simplified somewhat, but actually, the expressions 
given here are as compact as any alternate form. 

We are thus finally in a position to write the influence 
coefficients in the expressions 




(18) 



(19) 
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( 20 ) 
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and we obtain the evaluations 
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APPENDIX 3 



SHRINK ASSEMBLED MULTILAYER CYLINDRICAL 
VESSELS UNDER PRESSURIZATION 

(This appendix was contributed by Prof. J. E. Brock. 
His notation has been changed so as to conform with that 
appearing elsewhere herein.) 




% 



Pig. 1 Dimensions and notations for shrink assembled 
multilayer vessel. 



125 



Figure 1 indicates the configuration to be analyzed. 
The radii are prior to assembly and pressurization. The 
interferences 




(1) 



are small compared to either radius on the right. Accord- 
ingly in the analysis which follows we may use the 
approximation 




( 2 ) 



Also, it is clear that interface pressures are related by 



V. 



■fxtl 



C< 



- ¥ 
— if 






Co 



( 3 ) 



In particular 




the internally and externally applied pressures, respectively. 

By applying Lame's formulas (cf. Section IX hereof) to 
the i*'* 1 shell, we find 




(5) 
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(6) 



(7) 







r i 

do 



( 8 ) 



Assuming a state of plane strain with e = const, for 
all layers, fundamental equations of elasticity lead to 




(9) 



Note that this is constant throughout the 



i th shell. 




3»* £; 



(10) 



is the outward deformation at the outside (see Fig. 1) and 
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f* 



(id 



is the outward deformation at the inside. The constants are 




( 12 ) 



(13) 



( 1*0 



(15) 



By presuming that there are no gaps between shells under 
any conditions, geometrical continuity leads to the relation 






(16) 
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Substituting equations (10) and (11) into (16) gives 



\ a>-A{P eo ' = S* 



( 17 ) 



l - (h % ) °° °j M~l) 



Transposing known quantities to the right and assembling 
the results in matrix form, we have the system 



im ■ w 



+ 1 



V t L 



(18) 



where 
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(?,- B,.,) 


_/ Vi. 
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(C*,-\x) 
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^3 (Cj-'Bj) -A 2 

0 A (C--B,) 



is a triple band square matrix. 
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• CO 



( 19 ) 





•irt-2 






Is a column matrix of unknown interface pressures. 







( 20 ) 



is a column matrix of the radial interferences. 









( 21 ) 



is a column matrix containing differences between Poisson’s 
ratios, and 




_ V 



o , o, 



CO o 




( 22 ) 



is a column matrix involving exterior and interior fluid 
pressures . 
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By use of equation (9), we can find the total axial 
tensile force F as 



where 



and 



F.rr 



* = I 



r Ga +2,r < p 




( 23 ) 



iV , ,V.,,2 



<3= ir > &r-(cr ^ 
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( 25 ) 



- <( ( (K ) ) v^; o ) (e co ^ ( v^- ^ 1 



is a row matrix. 
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Equation (23) permits expressing e in terms of F as 

A 

follows 



F-fe-air <M>{P] 

G 



( 27 ) 



If e is given, there is no difficulty in solving 

A 

equation ( 18 ) and the solution may be written 






( 28 ) 



where 

fl-tsl 1 



(29) 



Having the interface pressures in {P} , it is elementary 
to find stresses o^ and in the individual layers. 

Equation (9) gives the axial stresses, and equation (23) 
gives the total axial force. 

If the net total tensile force F is given, equation (27) 
is used to replace e on the right side of equation ( 18 ), 

A 

which can then be written 






L 



(30) 
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where 



[S , l = ts],f{v] < M> 

is a square matrix which, however, is no longer banded, 

M - 1 4 

and 



The solution is 




i— n 



C A tf C* V 



•t 



C : 



1 * 



( 34 ) 



where 

[C-] - [sf 



( 35 ) 
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The interface pressures are now known and e may be 

A 

obtained from equation (27). Thus the stresses may easily be 
calculated. It should be remarked that this plane strain 
analysis gives results which are significantly different 
from results obtained by analyzing the assembly and loading 
of thin rings . The latter analysis is sometimes incorrectly 
given as being applicable to the case of multilayer cylinders. 
It should also be noted that the present analysis assumes 
that e from initial unpressurized unassembled condition 

A 

to final pressurized assembled condition, is the same for 
each cylinder. Practical methods of assembly usually imply 
adding one layer at a time and during the successive addition 
of layers to achieve the final assembled vessel, it may not 
be possible to assure that e is indeed the same for each 

A 

layer. However, after assembly is complete, the additional 
e due to pressurization is indeed the same for each layer. 

A 

The special case n = 2 is of particular interest and 
leads to the following evaluation if e is specified. 

A 








( e<'j J 






■ V -V‘ 

v \ y l 




[tey-ietf] J 



( 36 ) 



This may be compared with results given in Section L. 
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